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Abstract 

We study a manifestly gauge invariant set of Schwinger-Dyson equations to determine the non- 
perturbative dynamics of the gluon and ghost propagators in d = 3 Yang-Mills. The use of the 
well-known Schwinger mechanism, in the Landau gauge, leads to the dynamical generation of a 
mass for the gauge boson (gluon in d = 3), which, in turn, gives rise to an infrared finite gluon 
propagator and ghost dressing function. The propagators obtained from the numerical solution 
of these nonperturbative equations are in very good agreement with the results of SU{2) lattice 
simulations. 

PACS numbers: 12.38.Lg, 12.38.Aw, 12.38.Gc 
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I. INTRODUCTION 



QCD in three space-time dimensions (QCDg for short) has received increasing attention 
in recent years, not only because it is the infinite-temperature hmit of its four- dimensional 
counterpart (QCD4), but also because, at zero temperature, these two theories, despite 
a number of important differences, seem to share a variety of important nonpertubative 
features [l-28| . 

QCDg differs from QCD4 in several aspects. For example, the fact that QCD3 lives in an 
odd- dimensional space allows the appearance of phenomena that are not possible in even- 
dimensional spaces, such as the parity violating gauge-boson masses from a Chern-Simons 
term . In addition, unlike QCD4, there is no linearly rising potential for the quarks. 

Moreover, given that in = 3 the square of the coupling constant has dimensions of mass, 
QCD3 is super-renormalizable, having a trivial renormalization CToup. Finally, there are no 
finite-action classical solitons in QCD3 {i.e., no instantons) [see |9| for a brief review]. 

On the other hand, both theories confine, display area laws for Wilson loops in the fun- 
damental representation, and develop nonperturbative vacuum condensates, such as Ti (Qfj); 
in fact, in d = 3 one can actually prove [11] the existence of a TY{Qfj) condensate, associated 
with the minimum of the zero-momentum effective action, simply on the hypothesis that 
the full theory possesses a unique mass scale (that of the gauge coupling). In addition, and 
more importantly for the purposes of the present work, both theories appear to cure their 
infrared (IR) instabilities through the dynamical generation of a gauge boson (gluon) mass, 
usually refereed to also as "magnetic" mass, without affecting the local gauge invariance. 



which remains intact 



291 ]. The nonperturbative dynamics that gives rise to the generation 



of such a mass is rather complex, and can be ultimately traced back to a subtle realization 
of the Schwinger mechanism 3Ch37[. 

The gluon mass generation manifests itself at the level of the fundamental Green's func- 
tions of the theory in a very distinct way, giving rise to an IR behavior that would be difficult 
to explain otherwise. Specifically, in the Landau gauge, both in d = 3, 4, the gluon propaga- 
tor and the ghost dressing function reach a finite value in the deep IR. However, the gluon 
propagator of QCDg displays a local maximum at relatively low momenta, before reaching a 
finite value at g = 0. This characteristic behavior is qualitatively different to what happens 
in d = 4, where the gluon propagator is a monotonic function of the momentum in the entire 
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range between the IR and UV fixed points 



It should also be mentioned that a qualitatively similar situation emerges within the 



"refined" Gribov-Zwanziger formalism 39|, |40|, presented in 4l|. In this latter framework 



the gluon mass is obtained through the addition of appropriate condensates to the original 
Gribov-Zwanziger action. 

Even though several aspects of QCD3 have been studied in a variety of works, the recent 
theoretical developments associated with the pinch technique (PT), together with the high- 
quality lattice results produced, motivate the detailed study of the entire shape of the gluon 
and ghost propagators in d = 3. Specifically, given that the gluon mass generation is a 
purely nonperturbative effect, in the continuum it has to be addressed within the framework 
of the Schwinger-Dyson equations (SDE). These complicated dynamical equations are best 
studied in a gauge-invariant framework based on the pinch technique (P T) 29l . |42|-|45|| . and 
its profound correspondence with the background field method (BFM) [46']. As has been 
explained in detail in the recent literature 47|, |48|, this latter formalism allows for a gauge- 



invariant truncation of the SD series, in the sense that it preserves manifestly and at every 
step the transversality of the gluon self-energy. 

In the present work we study the dynamics of the gluon and ghost propagators of pure 
Yang-Mills in d = 3, using the SDEs of the PT-BFM formalism in the Landau gauge. Even 
though our results are valid for every gauge group, we will eventually focus on the group 
SU{2), in order to make contact with available lattice simulations The crucial ingredi- 
ent in this analysis, which accounts for the type of solutions obtained, is the gauge-invariant 
introduction of a gluon mass. The way gauge invariance is maintained is through the in- 
clusion of Nambu-Goldstone-like (composite) massless excitations into the non-perturbative 



three-gluon vertex 



29| . As a result, the fundamental Ward identities of the theory, which 



encode the underlying gauge symmetry, remain intact. The results obtained from our SDE 
analysis, presented in section 4, compare rather well with the available lattice data [see in 
particular Figs |6] and [8] . 

In addition, as a necessary intermediate step, we calculate an auxiliary function, denoted 
by G{q), which plays an instrumental role in the PT-BFM framework (see next section). 
Interestingly enough, and in the Landau gauge only, G{q) coincides with the so-called Kugo- 
Ojima (KO) function; this latter function, and in particular its value in the deep IR, is 
intimately connected to the corresponding and well-known confinement criterion 50 1. 
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The article is organized as follows. In Section II we briefly review the salient features of 
the SDEs within the PT-BFM framework. Section III contains a general discussion of the 
main conceptual issues related with the dynamical mass generation through the Schwinger 
mechanism. Particular attention is paid to the speciflc form of the three-gluon vertex that 
must be employed in order to maintain gauge invariance, in the form of the Ward identities. 
In addition, we give a qualitative discussion of some of the main features expected for the 
gluon propagator in the presence of a gluon mass. Section IV contains the main results 
of this work. After setting up the corresponding SDE for the gluon propagator and the 
auxiliary function G{q), we give explicit closed expressions for the latter quantities. The 
two available free parameters appearing in the expression for the gluon propagator, namely 
the gauge coupling g and the mass m are then varied, in order to obtain the best possible 
agreement with the lattice data. The ghost dressing function is also obtained from the 
self-consistent solution of the corresponding SDE; it too shows a good agreement with the 
lattice. Finally, in Section V we present our conclusions. 



II. THE PT-BFM FRAMEWORK 



In this section we remind the reader the basic characteristics of the SD framework that 



is based on the PT-BFM formalism; for an extended review of the subject see 44 1. 

We start by introducing the necessary notation. The gluon propagator ^^^{q) in the 
covariant gauges assumes the form 



(2.1) 



where ^ denotes the gauge-flxing parameter, Pfiv{q) = g^u — QfiQu/q"^ is the usual transverse 
projector, and A"^(g) = q'^ + ill{q), with n^j,y(g) = P^;y(g)n(g) the gluon self-energy. We also 
define the dimensionless vacuum-polarization n(g), as n(g) = g^n(g). In addition, the full 
ghost propagator, D{p), and its dressing function, F{p) , are related by D{p) = iF{p)/p^. 

The truncation scheme for the SDEs of Yang-Mills theories based on the PT respects 
gauge invariance {i.e., the transversality of the gluon self-energy) at every level of the 
"dressed-loop" expansion. This becomes possible due to the drastic modifications imple- 
mented to the building blocks of the SD series, i.e., the off-shell Green's functions themselves, 



following the general methodology of the PT 



29 



42. 



45| . The PT is a well-defined algorithm 
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FIG. 1: The full SDE for the gluon self-energy in the PT-BFM framework. By virtue of the special 
Abeliean-like Ward identities satisfied by the various fully dressed vertices, the contributions of 
each block are individually transverse. 

that exploits systematically the BRST symmetry in order to construct new Green's functions 
endowed with very special properties; in particular, the crucial property of gauge invariance, 
for they satisfy Abelian Ward identities instead of the usual Slavnov- Taylor identities The 
PT may be used to rearrange systematically the entire SD series 471]. In the case of the 
gluon self-energy it gives rise to a new SDE, shown schematically in Fig. [TJ 

Note that the quantity that appears on the Ihs of Fig. [T]is not the conventional self-energy 
H^u, but rather the PT-BFM self-energy, denoted by II^j,. The graphs appearing on the rhs 
contain the conventional self-energy II^i, as before, but are composed out of two types of 
vertices: 



(i) The conventional vertices, where all incoming fields are quantum fields, i.e., they carry 
the virtual loop momenta; these vertices are all "internal", i.e., the external gluons 
cannot be one of their legs, and will be generally denoted by T. 

(ii) A new set of vertices, with one of their legs being the external gluon, carrying physical 
momentum q; these new vertices, to be generally denoted by T, correspond precisely 
to the Feynman rules of the BFM {4^, i.e., it is as if the external gluon had been 
converted dynamically into a background gluon. 



As a result, the full vertices r;^™"(g, fci, fca), r^"™(g, fci, ^2 



' 5 a/iup 



{q,ki,k2,h), and 



T'^"'^ {q, ki, k2, ks) appearing on the rhs of the SDE shown in Fig. [T] satisfy the simple 




FIG. 2: Diagrammatic representation of the functions A and H. 



Ward identities 

g"f — = gr'^T^r^l + ^2, A:3) + cp., (2.2) 

where "cp" stands for "cyclic permutations" . Using these identities, it is straightforward to 
show that the crucial transversality condition q^Ii^y[q) = is enforced "block- wise" j4^, 



g^[(ai) + {a2%u = 0, 

g^Kas) + {a4)]^,u = 0, 

g^[(a5) + (ae)]/.. = 0, 

q^'iiar) + (og) + (ag) + {aw)]^.u = , (2.3) 

which allow for a self-consistent truncation of the full gluon SDE given in Fig. [1] 

Quite interestingly, the conventional A(g) and its PT-BFM counterpart A(g) the two 



5l| 



quantities are connected by the following background- quantum identity 

A(g) = [l + G(g)f A(g), (2.4) 
where the function G{q) is the g^i, component of the auxiliary two-point function A^^{q), 
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q q q k+q 

FIG. 3: The SDE satisfied by the ghost propagator. 



defined as 



Ca f Hf^D{k + q)^P"{k)HUk.q), 

Jk 



= 9,Mq) + ^-f^m, (2.5) 

where Ca the Casimir eigenvalue of the adjoint representation [Ca = for SU{N)] and 
Jj^ = fi'^'^{27T)^'^ Jd'^k, with d the dimension of space-time. The function Hcru is given 
diagrammatically in Fig. O Note that it is related to the full gluon-ghost vertex by 
q''H„i,{p,r,q) = -iT^{p,r,q)\ at tree-level, h'^] = ig^^. 

The identity (12. 4p allows to express the SDE of Fig. [1] as an integral equation involving 
only A(g), namely 



^-\q)PAq) = ^'-^^-y+;^^;=/^-^^- . (2.6) 



^10 

[1 + G{qW 

Finally, as shown in Fig. |3l the ghost SDE is the same as in the conventional formulation, 
namely 



tD'\q) = q^ + iq^Ca [ T^A^,{k)T''{q, k)D{q + k), (2.7) 

Jk 

where F^ is the standard (asymmetric) gluon-ghost vertex at tree-level, and it fully- 
dressed counterpart. 



III. MASS GENERATION IN d = 3 YANG-MILLS 

It is well-known that, just as happens at = 4, the Yang-Mills dynamics in c? = 3 
generates an effective gauge-boson mass, which cures all IR instabilities. The underlying 
mechanism that leads to the generation of such a dynamical mass, both in d = 3, 4, is the 
Schwinger mechanism, the only known procedure for obtaining massive gauge bosons while 
maintaining the gauge-symmetry intact. 
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FIG. 4: Vertex with non-perturbative massless excitations triggering the Schwinger mechanism. 



As Schwinger pointed out long time ago js^, the gauge invariance of a vector field does 
not necessarily imply zero mass for the associated particle, if the current vector coupling 
is sufficiently strong. According to Schwinger 's fundamental observation, if 11 (g) acquires 
a pole at zero momentum transfer, then the vector meson becomes massive, even if the 
gauge symmetry forbids a mass at the level of the fundamental Lagrangian. Indeed, it is 
clear that if the vacuum polarization 11 (g) has a pole at = with positive residue /x^, 
i.e., n(g) = /i^/g^, then (in Euchdean space) A~^(g) = g^ + fi'^. Thus, the vector meson 
becomes massive, A^^(O) = /i^, even though it is massless in the absence of interactions 
{g = 0). There is no physical principle which would preclude n(g) from acquiring such 
a pole, even in the absence of elementary scalar fields. In a strongly- coupled theory, like 
non-perturbative Yang-Mills in d = 3,4, this may happen for purely dynamical reasons, 
since strong binding may generate zero- mass bound-state excitations 32j. The latter act 
like dynamical Nambu-Goldstone bosons, in the sense that they are massless, composite, 
and longitudinally coupled; but, at the same time, they differ from Nambu-Goldstone bosons 
as far as their origin is concerned: they do not originate from the spontaneous breaking of 



any global symmetry 



29j. In what follows we will assume that theory can indeed generate 



the required bound-state poles; the demonstration of the existence of a bound state, and 
in particular of a zero-mass bound state, is a difficult dynamical problem, usually studied 
by means of integral equations known as Bethe-Salpeter equations (see, e.g.,{52|). Note also 
that the generation of a dynamical mass (both in d = 3,4) requires (and, correspondingly, 
gives rise to), the formation of a gluon condensate. 

The Schwinger mechanism is incorporated into the SDE of the gluon propagator essen- 
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tially through the form of the fully-dressed, nonperturbative three-gluon vertex. In fact, 
since the generation of the mass does not interfere with the gauge symmetry, which remains 
intact, the three-gluon vertex must satisfy the same Ward identity as in the massless case 
[viz. Eq. fl2.2p ]. but now with massive, as opposed to massless, gluon propagators on its 
rhs. The way this crucial requirement is enforced is precisely through the incorporation into 
the three-gluon vertex of the Nambu-Goldstone (composite) massless excitations mentioned 
above. To see how this works with a simple example, let us consider the standard tree-level 
vertex 

^tiai3{q,P, r) = {q- p)pg^,a + (p - r)^^„/3 + (r - q)a9^,f} , (3.1) 
which satisfies the simple Ward identity 

q^T,^p{q,p,r) = P^p{r)A,\r) - P^^{p)A,\p) (3.2) 

where AQ^(g) = is the inverse of the tree- level propagator. After the dynamical mass 
generation, the inverse gluon propagator becomes, roughly speaking, 

A;;,\q) = q'-m'iq'), (3.3) 

and the new vertex, r™^^(g,p, r) that replaces T^ai3{q,P, r) must still satisfy the Ward iden- 
tity of (13.21) . but with Aq ^ — on the rhs. This is accomplished if 

^^^a^iQ^P^r) = r^a/3(g,p,r) + V^«^(g,p,r), (3.4) 

where Vij^ap{q,p,r) contains the massless poles. A standard Ansatz for V^a/3(5',p, t) is [3] 

V,apiq,p, r) = m2(r)^^^^||^ P;(r) - [m^p) - m'iq)] ^ P;(g) P^ip) + c.p. , 

It is easy to check that 

q^V^cxi}{q,p,r) = Pap{p)w?{p) - Pap{r)w?{r) , (3.5) 

and cyclic permutations. Therefore, one has 

g'^r- ^(g,p,r) = P^p{r)A-^\r) - P^p{p)A-^\p) , (3.6) 

as announced. Note that for constant masses [m(g) = m{p) = m{r) = m] the vertex of (13. 5p 
reduces to 



V^a(siq,P,r) = — 



\qi^pM-p)p ^pf^, , Parp{p-r)p ^^^^^ ^ r^qf,{r-q)p ^^^^^ ^^^^ 



-^-T— + P;(g) + ^^^^^^ PM 



Even though the precise implementation at the level of the complicated integral equations 
is rather subtle, the final upshot of introducing a vertex such as F™ (or more sophisticated 
versions of it) into the SDE for the gluon self-energy is that one finally obtains, gauge- 
invariantly, a non- vanishing A^^(O) and A^^(O). Qualitatively speaking, in Euclidean space 
and d space-time dimensions, the (background) gluon propagator is given by 

A-i(g) = g2 + n(g) + A-i(0), (3.8) 

where n(g) has the general form 

n(g) = ci/ / A(fc)A(A; + q)K^{q, k) + c^g^ [ D{k)D{k + q)K2{q, k). (3.9) 

Jk Jk 

The functions Ki{q, k) and K2{q, k) are SD kernels, whose closed form depends, among other 
things, on the dimensionality of space-time, the details of the vertices employed, and the 
gauge chosen, as do, in general, the constants ci and C2. Setting A^^(O) = m^, one then 
obtains 

A-^(g) =g2 + m2 + n(g). (3.10) 

To obtain the perturbative (one- loop) expression for A(g) one must substitute in the integral 
on the rhs of f l3.9p the tree-level values for A, Z), Ki and K2-, which is a good approximation 
for large values of the physical momentum q. However, for low values of g, one must solve 
the integral equation, which, under suitable assumptions, will furnish massive (IR finite) 
solutions for A(g). 

An easy way to qualitatively appreciate the effect of the mass on the solutions for A(g) 
is to substitute A — t- A^ in the first integral on the rhs of fl3.9p . assuming for simplicity 
a constant mass m, and use tree-level expressions for all other terms. This will furnish an 
approximate expression for n(g), to be denoted by 11^(5'), and the resulting A^^(g) will 
read 

K-\q) = q^ + m' + IL^{q). (3.11) 
In d = 4 the corresponding Iim[q) will have the form 

^miQ) = bg^Q^ [ dx\n[q^x{l~x)+m^]. (3.12) 
Jo 

For m — 0, or ^ m^, one recovers the usual one-loop logarithm bg'^ln{q'^), with b being 
the first coefficient of the QCD one-loop /3-function, b = 11Ca/4:8it'^. As explained in the 



10 



literature, the presence of the mass inside the logarithm tames the Landau pole, and gives 
eventually rise to a IR finite value for the QCD effective charge 
Similarly, in li = 3 we have [see the integral Ri in Eq. fl4.8p ] 

&^\q) = -2hg'qaTctan (^) , (3.13) 

which in the limit m — )■ assumes the one loop perturbative form [see also the integral Ii 
in Eq. 1^' 

n£lt(g) = -^h9\ (3.14) 

In this case however, and unlike in ci = 4, 63 is a numerical coefficient that depends explicitly 
on the value of the gauge parameter chosen; in the Feynman gauge, 63 = 15CA/327r (we will 
return to this point in the next section). 

Let us now briefiy compare the versions of the gluon propagator obtained by substituting 
npprt(g) or Urn (q) into (13. lip . For the perturbative case we have 

Apert(g) = 2 \ 2 ■ (3-15) 

q - T^hg q 

There two points to notice: (i) Apert(5') has a Landau pole at g = nb^g'^, and (ii) it displays 
a maximum value at q* = q/2. On the other hand, the gluon propagator corresponding to 
^m\q) becomes 

It is clear that the presence of the mass regulates the denominator for all values of q, provided 
that it exceeds a certain critical value (in units of g'^). In addition, A(g) may or may not 
display a maximum, depending on the ratio g'^/m; in general, its position is displaced with 
respect to q*. 



IV. RESULTS AND COMPARISON WITH THE LATTICE. 



In order to make contact with the d = 3 lattice results of 49||, we must next determine 
the form of the relevant SDEs in the Landau gauge (^ = 0). The three quantities of interest 



are: 



(i) The gluon propagator, A(g) given in (12.60 : 
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(ii) The Kugo-Ojima function G{q), given in (12. Sp . which connects the conventional and 
background gluon propagators; 



(iii) The ghost propagator, given in (12. 7p . and in particular its dressing function, F{q). 



A. Calculating the gluon propagator(s) and the KO function 

In the "one-loop dressed" approximation, the PT-BFM gluon self-energy is given by the 
following (gauge- invariant) subset of diagrams: 

= + («2) + (as) + (04)]^^. (4.1) 



When evaluating the diagrams (a^) one should use the BFM Feynman rules 46|, noticing in 
particular that the bare three- and four-gluon vertices depend explicitly on 1/^, the coupling 
of the ghost to a background gluon is symmetric in the ghost momenta, and finally that there 
is a four-field coupling between two background gluons and two ghosts. 

As explained in [38|, the limit ^ — > of the diagrams (oi) and (02) must be taken with 
care, due to the terms proportional to 1/^ coming from the tree-level vertices. Introducing 
^Uil) = Pf^u{q)Mq), one obtains 

[{a,) + {a,)U = /C'A|iy^r^'^A*,^(A:)A^^(fc + g)C-^^7,.y^A(A:) 

The vertex l^ai3 is the fully-dressed counterpart of T^a/s (in the Landau gauge); it satisfies 
the Ward identity 

g%a/3 = Pai^ik + q)A-\k + q)- P^^{k)A-\k). (4.3) 

It is then easy to verify that the rhs of (14. 2 p vanishes when contracted by q'^, thus explicitly 
confirming the validity of the first equation in ( 12. 3p . for the special case of ^ = 0. 
Similarly, 



[(as) + {ai)]^u = -g^C/ 



I T^D{k)D{k + q)T, - 2ig^, I D{k) 

Jk Jk 



(4.4) 



with r^(g,p, r) ~ (r ~p)p- The vertex satisfies the second Ward identity in (12. 2p . which 
leads immediately to the transversality of this block, i.e., the second equation in (12. Sp . 
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Finally, using tree-level values for the auxiliary function H^jy in (12.50 and for the vertex 
in (12 .yp . we obtain for the Kugo-Ojima function 



G{q) 



while for L[q) one has 



L{q) 



1 + 



1-3 



A{k)D{k + q), 



A{k)D{k + q). 



(4.5) 



(4.6) 



The way we proceed is the following. Instead of actually solving the system of cou- 
pled integral equation, we will adopt an approximate procedure, which is operationally less 
complicated, and seems to capture rather well the underlying dynamics. 

Specifically, we will assume that the gluon propagator has the form given in (13. lip , and 
will determine the function Hm (q) by calculating the expressions given in ( 14. 2p and (14. 4p 
using inside the corresponding integrals A — )■ and D ^ Dq. In order to maintain gauge 
invariance intact, we will set 



= r;r,^(?,p,0 (4.7) 

with rj^^(g,p, r) given in (13. 4p . The vertex V^^(g,p, r) entering into r^^(g,p, r) will be 
that of Eq. (13. 7p . i.e., we will assume a constant mass m throughout. 

From the final expressions appearing in the rest of the paper we will use Euclidean 



momenta. To that end we set 



with g| > the positive square of a Euclidean 



four-vector, and = ygl- The Euchdean propagator is defined as A£;(g|) = — A(— g|). 
To avoid notational clutter, we will suppress the subscript "E" in what follows. 

The results of all our calculations will be expressed in terms of the following six basic 
integrals, 

1 (± 

1 f i \ I 



Ro 
Ri 
h 
h 
h 
h 



f^2 _ ^2 



(A;^ — m^)[(A; + g)^ — mP\ 
1 ( i\ 1 

sj g' 

i 



A.TX ) q 



arctan 



q 

2m 



kk'ik + qy 
1 



k (k^ - m^)ik + q) 

1 f i \ I 

An ) m 



- arctan {— 
An J q \m 



^ fc2(A;2 — m^) 
g ■ k 



(F - m2)(A; + qy 



m + 



m 



arctan ( — 



(4.8) 
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where the momentum q appearing in the integrals on the Ihs is Minkowskian, while the 
momentum q appearing in the results on the rhs is Euclidean. 

To facilitate the calculation, and since the transversality of Ii^y{q) is guaranteed, one 
may set in (14. ip Ii^y{q) = P^^{q)Ilm{q), and isolate Umil) by taking the trace of both sides, 
i.e., 

{d - l)flm{q) = [(ai) + (as) + (as) + (a^)]'; (4.9) 
For the different four contributions shown in Eq. (14. 2 p we obtain the following results 

I [ TfAi^ik)Al^{k + q)F^'^ = 9Ro + f ^4 -^—2- 10?' + 8^') + 
I Ji^ ^ ^ ^ \4m^ m"^ / 4m^ 



2 m? 2 ^ / V 2 



'3, 



k 



3) 



A;nA; + g)^ 1 



2r 

- 'g'/i. (4.10) 



+ ( 1^ _ 2ig2 _ 3^2 + ±11^ J + ± 2 ^ ::_ I J. 
\2m2 4 ^ 4 g2 y ^ \^2^ 4 ' ' 



lkk^{k + qy 2 

Next, let us turn to the diagrams (as) and (a4) of Fig. [H which contain a ghost loop. Since 
we will treat the ghost as a massless particle, the "tadpole" diagram (a4) vanishes identically 
in dimensional regularization; from diagram (as) we get instead (after taking the trace) 



("3)2 - -ic. I ^ -/c.,=A. (4.11) 



k 



From the results above it is relatively straightforward to check, taking appropriate limits, 
that in the deep IR 

n(0) = -il-^m. (4.12) 
Qn 

Therefore, in order for the (Euclidean) A(0)^^ = — in(0) to be positive definite, m and 
g must satisfy the condition 

ml , , 

> — . (4.13) 



In the opposite limit, namely for asymptotically large momenta, the addition of all terms 
given in fl4.10p exposes a vast cancellation of all powers g", with n > 1. After all such 
cancellations taking place, one is left with a linear contribution, given by 

n(g) f 15 - ly (4.14) 



32 n 2 
14 



The reason for writing the numerical coefficient in front of the leading contribution as a 
deviation from 15 is the following. The expression (I4.14p should coincide with the d = 3 
one-loop BFM self-energy calculated in the Landau gauge. For any dimension d and any 
value of the gauge-fixing parameter C,q, the latter reads 44 1 



d-1 Jk k^{k + qY 



^ "^^^ ]uk\k + qY^ J, 



k\k + qY 



(4.15) 



In the Feynman gauge of the BFM, = 1, n(g) collapses to the PT answer for the gauge- 
independent gluon self-energy; specifically, for d = 3, 

n(g)k,=i = -^^?(15). (4.16) 

Away from = 1 the terms in the second line of f l4.15p give additional contributions, which 
may be easily calculated using the basic results 

kp^ky i 1 

kk^ikTW " ~32'^^^ ' 

= - J- (4 17) 

Ikk^k + qf 16g' ^ ' ' 

where the dots in the ffist integral indicate longitudinal parts. In particular, it is easy to 
verify that at = these additional terms account precisely for the term — | appearing in 
Eq. dUl. 

The above discussion reveals an important difference between the d = 3 and d = 4 cases. 
Specifically, in ci = 4 the coefficient in front of the leading one-loop contribution to n(g) is 
independent of the gauge- fixing parameter ^g. This well-known BFM result can be easily 
deduced from fl4.15p . since both integrals proportional to (1 — ^g) are UV finite, i.e., they do 
not furnish logarithms. The coefficient in front of the logarithm is completely determined 
by the ffist integral, multiplied by the factor (itt)? which, at (i = 4, reduces to 
1Qtx% = (11/3)(7^Ca, namely the ffist coefficient of the Yang- Mills /3 function. As we have 
just demonstrated, things are different in c? = 3, where no renormalization is needed; the 
leading (linear) contribution depends explicitly on the value of ^q. 

Next, we determine an approximate expression for the function G{q). To that end, we 
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turn to (14. 5 p and substitute in the integral on the rhs, A — t- Am and D ^ Dq. One has then 



G{q) 



1 + 



k {k + qy{k'^ — m?) [ ' k'^q^ 

■-/4 + 5/2-/3 + ^(/2-/l) 



(4.18) 



which gives 



327rm 



TT q m" I m? q^ 

H 7^ — 1+6 I arctan 



2 m q^ 



(4.19) 



In the deep IR (g — ?■ 0), and for asymptotically large momenta {q — t- oo), one finds 



G{q) 



G{q) 



0. 



(4.20) 



From the expressions for Tlrniq) and G{q) obtained above, we can determine the conven- 
tional gluon propagator, A(g), in the Landau gauge; the latter can then be compared to the 
lattice data. To that end, let us first employ the crucial identity of (12. 4p to write 

[^ + G{q)f 



A{q) = 

q^ + 1V? + Iim{q) 

Then, by virtue of (I4.20p . A(g) has the same asymptotic behavior as A(g). 



(4.21) 



Notice that in d = 3 the gluon and ghost propagators have the basic scaling property 



A(g, g, m) = a^A{aq, y/ag, 



am) 



D{q, g, m) = a^D{aq, \fag^ am), 



(4.22) 



where a is a positive real number. Of course, the corresponding dressing functions (being 
dimensionless quantities) are invariant under such a combined rescaling; for example, the 
ghost dressing function satisfies F{q, g, m) = F{aq, y/ag, am), and so does the gluon dressing 
function g^A(g) and the Kugo-Ojima function G{q). One can then make use of these scaling 
properties to set g (respectively m) equal to unity, and vary m (respectively g) in order to 
study the shape of the solutions found so far. The results (when setting g = I and varying 
m) are shown in Fig. [5l 



B. Comparing the gluon propagator with SU{2) lattice data 



We next compare the result of our calculation for the conventional gluon propagator 
A(g) with the lattice results of j49|. In order to do that, the lattice data must be first 
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FIG. 5: Results for the massive one-loop approximation for the d = 3 gluon propagator. In 
the upper panels we show the plots for different values of the hard-mass parameter m for the 
background-quantity identity ingredients A(g) (left) and the Kugo-Ojima function G{q) (right). In 
the lower panels we show the conventional propagator A{q) (left) and its corresponding dressing 
function q'^A(q) (right). 

properly normalized (or, equivalently, the theoretical prediction must be suitably rescaled) 
Specifically, in the absence of any physical input that would fix the physical scale, one 
uses the scaling property f l4.22p and determine the scaling factor a in such a way that the 
asymptotic (large momentum) segment of the lattice data coincides with that obtained from 
our calculation; indeed the two "tails" should coincide, given that perturbation theory is 
reliable in that region of momenta. The result of this procedure is shown in Fig. [6l evidently, 
the matching between the theoretical curve and the lattice data is very good. The best-fit 



17 



0.8 
0.6 
0.4 
0.2 






P=2.2 ^Aol 




200, — a - 








320-' 




Theory 







0.5 1 1.5 

q [lattice units] 




0.5 1 1.5 

q [lattice units] 



FIG. 6: Comparison of the lattice results of [49] with the gluon propagator (left) and the gluon 
dressing function (right) obtained within the massive one-loop approximation adopted in this paper. 
In passing, notice that the dressing function does not tend to 1 for asymptotically large q which 
also motivates the momentum rescaling procedure employed. 



curve furnishes the ratio 



m 
2? 



0.146, 



(4.23) 



which appears to be in rather good agreement with previous theoretical and lattice stud- 



ies 



y. 



C. Ghost dressing function and lattice data 

We next proceed to calculate the theoretical prediction for the ghost-dressing function 
F{q). In a spirit similar to that adopted for the gluon propagator, as first approach in this 
direction, we simply compute the diagram for the ghost propagator (see Fig. [3]) using as 
inputs on the rhs A — )■ A^^ and D ^ Dq. The result of this calculation is 



F{q) 



1 + 
1 + 



1 - 



J. {k + qY{k'^ — rn?) 

+ 3/2 + /s - ^ {h 



k'^q^ 
h) 
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(4.24) 



At this point, and before attempting a comparison with the corresponding lattice data, we 
note that, in the Landau gauge only, the ghost-dressing function F{q), and the two form 
factors G{q) and L{q) defined in Eq. (12.51) . are related (for any d) by the following important 
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identity, 



l + G{q) + L{q) = F~\q). 



(4.25) 



The relation of Eq. fl4.25p . has been first obtained in 53|], and some years later in in 
the framework of the Batalin-Vilkovisky quantization formalism; as was shown there, this 
relation is a direct consequence of the fundamental BRST symmetry. Recently, the same 



identity has been derived exactly from the SDEs of the theory [55j, and the important 
property L{Qi) = 0, usually assumed in the literature, was shown to be valid for any value of 
the space-time dimension d; indeed setting A — )■ Am and D ^ Dq on the rhs of Eq. (14.61) . 
one has 



Hq) 



^9^Ca 
327rm 



1 



k{k + qYik?' — m?) 



TT q m 

-— H 

2 m 



1 + 



m 



m^ 



arctan 



m 



and therefore 



0, 



L{q) 0. 



(4.26) 



(4.27) 



It is then straightforward to verify from the result above and the closed expressions given 
in Eqs (14.191) and (14.241) . that Eq. (14.251) holds exactly within the approximation scheme we 
are using (see also the left panel of Fig. [7]). 

We next vary the parameters g and m in the expression given in Eq. (I4.24p in order to 
reproduce the lattice data for F{q). As a natural starting point we use the values that have 
resulted in the best fit for the gluon propagator, namely g = 1.285 and m = 0.480. However, 
as is clear from the red dashed curve shown in Fig. [7] (right panel), the result obtained is in 
poor agreement with the lattice. If instead we allow the parameters to vary freely, i.e., we 
disregard the gluon data and attempt to only reproduce the ghost data, the best possible 
curve is shown by the black continuous line of Fig. U\ being obtained for the values g = 2.049, 
m = 0.543 (giving m/2g'^ ^ 0.065). 

It is clear from this analysis that, within the approximation scheme employed, the lattice 
data may be well reproduced if treated independently, but it is not possible to arrive at a 
reasonable simultaneous fit, i.e., to fit both curves using a unique set of parameters. 
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FIG. 7: (Left panel) Values of L{q), 1 + G{q), and F{q) = [1 + G(q) + L{q)\^'^ within our approx- 
imation for g = 1.287 and m = 0.539. (Right panel) Comparison of the ghost dressing function 
with the one calculated within our approximation for two sets of values corresponding to the gluon 
propagator best fit {g = 1.287 and m = 0.539, red dashed line) and to the best fit to the lattice 
data {g = 2.049, m = 0.543). 

D. Combined treatment: gluon propagator and ghost dressing function vs lattice 

To remedy this situation, we will improve the approximation used for obtaining the theo- 
retical prediction for the ghost-dressing function. Specifically, we will study an approximate 
version of the ghost SDE given in Eq. (12. 7p . and we will solve self-consistently for the un- 
known function -F(g), instead of simply calculating its rhs, as was done above for obtaining 
the expression in Eq. (14. 24 p . 

Given that Eq. (12. 7p contains ^{k) as one of its basic ingredient, the general matching 
procedure becomes more subtle. In particular, instead of freely fitting just one set of data 
(that for the gluon propagator) one must now attempt to fit simultaneously both the gluon 
and ghost data, as well as possible. As we will see, this more complicated procedure furnishes 
finally a very good agreement with the combined set of lattice data, but one has to settle 
for a slightly less accurate description of the gluon data compared to the one obtained in 
Fig.E] 

After approximating the gluon-ghost vertex T ^ by its tree-level value, we arrive at the 
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FIG. 8: Comparison of the lattice results of 49|] with the ghost dressing function obtained through 
the solution of the ghost SDE (left); on the right we show the gluon propagator for the same value 
of m/2c/2 0.153. 



following integral equation for the ghost dressing function F, 

(k-q)'^] A{k)F{k + q) 



F-\q) = l + g'CA 



{k + qf 



(4.28) 



The general idea now is to solve Eq. fl4.28p numerically for F{q), using as input for the A[k) 
under the integral sign the theoretical curve that, after the rescaling mentioned earlier, 
provides the best possible fit to the gluon data, and, at the same time, allows for the 
numerical convergence of Eq. fl4.28p . We note in passing that this procedure permits, after 
a shift of the integration variable, to pass all angular dependence from F{k + q) to A(fc), 
whose functional form is considered known; as a result, one does not need to resort to further 
approximations for the angular part of the integral equation. 

The general observation regarding the numerical treatment of Eq. fl4.28p is that it appears 
to be extremely sensitive to the precise shape of A(A;) and the value of g; minute variations 
of these quantities give rise to large disparities in the resulting F{q). 

The best possible solution that we have obtained is shown in the left panel of Fig. M As 
announced, the accuracy achieved in matching the lattice data for the gluon propagator is 
slightly inferior to that of our best fit (Fig. [HI right panel), but is still very good. 
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V. CONCLUSIONS 



In this work we have presented a nonperturbative study of the (Landau gauge) gluon and 
ghost propagator for d = 3 Yang-Mills, using the "one-loop dressed" SDEs of the PT-BFM 
formalism. One of the most powerful features of this framework is that the transversality 
of the truncated gluon self-energy is guaranteed, by virtue of the QED-like Ward identities 
satisfied by the fully-dressed vertices entering into the dynamical equations. 

The central dynamical ingredient of our analysis is the assumption that the famous 
Schwinger mechanism, namely the dynamical formation of zero-mass Nambu-Golstone- 
boson-like composite excitations, which allow the gauge-invariant generation of a gauge- 
boson mass, is indeed realized in d = 3 Yang-Mills. The way this dynamical scenario is 
incorporated into the SDEs is through the form of the three-gluon vertex. Specifically, in 
order to satisfy the correct Ward identity, as required by gauge- invariance, this vertex must 
contain massless, longitudinally coupled poles, representative precisely of the aforementioned 
composite excitations. 

It should be emphasized that the approach followed here is approximate, not only in the 
sense that we consider the one-loop dressed version of the SDE, omitting (gauge-invariantly) 
higher orders [i. e., the third and fourth block of Fig. [1], but also because we do not actually 
solve simultaneously the full system of resulting equations. Specifically, as explained in 
Section IV, when evaluating the gluon self-energy we have used tree-level expressions for the 
ghost propagators appearing in diagram (03) of Fig. [1], and the same approximation is used 
also in the determination of the KO function G{q). We have then used the resulting gluon 
propagator as an input into Eq. fl4.28p to obtain the improved F{q). Of course, this two-step 
procedure is bound to result in a considerable discrepancy between the "one-loop" G{q) and 
the F{q) obtained from solving its corresponding SDE; evidently, the identity of Eq. f l4.25p 
cannot be fulfilled any longer. In addition, the dynamical gluon mass m has been treated 
for simplicity as a constant. However, a more thorough study should eventually include the 



important feature that 
general considerations 



;he mass depends nontrivially on the momentum, in accordance with 



56| : in fact, a complete SDE treatment ought to actually determine 



the precise way the mass is running |57|. The fact that, despite these simplifications, the 
lattice results for the gluon and ghost propagator are so well reproduced, suggests that the 
full treatment may reveal a number of subtle cancellations, caused by the highly non-linear 
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nature of the SDE equations, yielding finally results very similar to those reported here. 

Let us finally outline briefly some of the modifications and additional field-theoretic inputs 
that such a full SDE treatment would entail. To begin with, a more complete Ansatz for the 
three-gluon vertex r^^u, appearing in the SDE of the gluon propagator [graph (ai) in Fig. [Tj, 
must be devised. Such an Ansatz must not contain only the part of the massless poles, as 
Eq. (13. 5p does, which only accounts for the massive part of the propagator, but should make 
explicit reference to the entire A, in the spirit of the analysis already presented in 57|] (for 
d = 4). In addition, a similar Ansatz must be introduced also for the full gluon-ghost vertex 
Fq, appearing in the graph (03) of Fig. [TJ In order to maintain explicit gauge- invariance. 
Fa must be such that the second Ward identity of Eq. (12.21) is automatically satisfied. Note 
that Ta is not the same as the conventional gluon-ghost vertex that appears in Eq. (12. 7p 
and in Fig. [31 Given that F^ and F^ are different, and that only the former is crucial for the 
transversality of the gluon SDE, one may approximate F„ by its tree-level value, without 
clashing with gauge-invariance, a freedom that exists only within the PT-BFM scheme. 
Should one opt for a more sophisticated treatment of F^, then an appropriate Ansatz may 
have to be devised. Given that F^ satisfies a complicated Slavnov- Taylor identity, instead 
of the simple Ward identity of F,^ , further approximations may be necessary. We hope to 
be able to implement some of the aforementioned improvements in the near future. 
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